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Bound-free e+e pair creation with a linearly polarized laser field and a nuclear field 

C. Deneke and C. Miiller 
Max- Planck- Institut fur Kernphysik, Saupfercheckweg 1, 69117 Heidelberg, Germany 

The process of bound-free pair production of electrons and positrons in combined laser and 
Coulomb fields is investigated. It is assumed that an ion collides at relativistic speed with an 
intense x-ray laser beam of linear polarization. The process proceeds nonlinearly due to simulta- 
neous absorption of a few laser photons. The capture of the electron into the ground state and 
the L-shell is considered. The scaling of the total rate, the angular distributions of the emitted 
positrons and a comparison to the competing free-free channel are surveyed. Numerical results of 
pair production rates for parameters for the planned x-ray free electron lasers at DESY and SLAC 
are presented. We find that pair production with these laser facilities can become observable in the 
near future. 

PACS numbers: 12.20.Ds, 32.80.wr, 34.90.-I-9, 42.55.Vc 
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I. INTRODUCTION 

Theoretical and experimental study of electron- 
positron pair creation in strong external fields is impor- 
tant for understanding the structure of the QED vacuum. 
In the past, pair production in relativistic heavy ion colli- 
sions was studied in detail 1, 2, 3]. With the progress of 
laser technology, pair production via absorption of sev- 
eral real photons becomes increasingly interesting. There 
are various schemes of pair production in laser fields, for 
example theproduction with counter-propagating laser 
fields (e.g. [J, [a, [g, 0, §1 and references therein). How- 
ever, for the observation it is necessary to come close to 
the critical electrical field strength £c = 1.3 x lO^^V/cm 
[9| which is still far away from experimental realisation. 

Another way is to consider pair production in com- 
bined laser and Coulomb fields. The Coulomb field can 
originate from a moving ion. If the ion countcrpropa- 
gates the laser beam at high speed, the laser's electric 
field strength and photon frequency are Doppler shifted 
by the factor « 27, where 7 is the relativistic Lorentz 
factor of the ion. One usually distinguishes between the 
multiphoton and the tunneling regime of pair produc- 
tion. In the latter, the laser field strength is close to the 
critical one and pair creation proceeds via the simultane- 
ous absorption of a very large number of low frequency 
photons (^ 10^). An important quantity which arises in 
strong- field physics is the intensity parameter, ^ = -£^^ 
with the laser frequency to and peak field strength Sq- 
Moreover, e, m and c denote the elementary charge, elec- 
tron mass and speed of light, respectively. The tunneling 
regime corresponds to ^ ^ 1, thus to a high intensity and 
low frequency laser field. Contrary, in the multiphoton 
regime the field strength is smaller, but the photon en- 
ergies are sufficiently high to allow pair creation via the 
absorption of only a few photons. Therefore, for the mul- 
tiphoton regime 1^ ^ 1. When electron-positron pairs are 
created in the vicinity of a nucleus, the electrons can be 
created either free or in a bound state. One refers usually 
to free-free and bound-free pair production, respectively. 
In the present work we consider bound-free pair creation 
in the multiphoton regime due to the head-on collision of 



a very fast bare ion or proton with a very intense X-ray 
laser field of linear polarization. 

X-ray free electron lasers (XFEL) currently under con- 
struction at DESY (Hamburg, Germany) and SLAC 
(Standford, California) are predicted to obtain photon 
energies of fiw ~ 4—12 keV and intensities of / '^ 
lO^^W/cm^ [1^. Energy conservation demands nhio > 
2mc^, the energy of n absorbed laser photons has to be 
greater than the rest masses of the produced particles. 
Hence, together with acceleration of ions to a relativis- 
tic 7- factor of around 50, pair creation with two or three 
photons from an XFEL is possible as a nonlinear process. 

In 1997, the process of laser-induced pair creation has 
been observed at SLAC [ill, |l2, Il3| • A highly relativistic 
electron beam collided with a highly intense laser beam. 
Energy conservation allowed the production of electron- 
positron pairs by absorption of five photons. Two differ- 
ent mechanisms were found to be responsible for the pro- 
cess. On the one hand the Breit- Wheeler process where 
in a first step one or two photons are absorbed to cre- 
ate a highly energetic photon via Compton scattering. 
In a second step this 7 photon interacts with additional 
laser photons to create pairs. Alternatively, pair produc- 
tion happens via the nonlinear Bethe-Heitler process in 
a single step by simultaneous absorption of several real 
photons in the Coulomb field of the electrons to create 
pairs. In the analysis of the experiment it was found that 
the first process dominated over the second. In contrast, 
when the laser beam collides with heavy ions instead 
of electrons, the high masses of the ions suppress the 
Breit- Wheeler process substantially and the (nonlinear) 
Bethe-Heitler process plays the dominant role. Bound- 
free pair creation has not yet been observed in laser-ion 
collisions, but in a similar way bound electrons have been 
detected in an experiment with relativistic heavy-ion col- 
lisions [iJi \lM, [3 ■ For small nuclear charge numbers Z 
bound-free pair creation is suppressed, but it is competi- 
tive with the free-free channel at intermediate and high Z 
Q . The two competing channels of free- free and bound- 
free pair creation could be distinguished experimentally 

Mm- 

Theoretically the process of free pair creation by com- 



bined Coulomb and laser fields has been investigated by 
several authors [H [S [H, [H [H [H, [H, [13, [M HI, M, 
[23, [23, [so, [3l| with a focus on the tunneling regime. Also 
muon pair creation in laser-ion collisions has been calcu- 
lated recently [32, ^] . Bound- free electron-positron pair 
creation with single high-en ergy photons has been sur- 
veyed before (see, e.g., [3J, [3^|), but for the nonlinear 
case, i.e. the absorption of more than one photon, only 
a few calculations exist [3a [33| . Both consider K-shell 
capture of the electron and a circularly polarized laser 
beam [3dl or an ultrashort, single-cycle electromagnetic 
pulse [37|, respectively. 

It is worth mentioning that in relativistic laser-ion or 
laser-electron collisions also other QED processes can 
take place which arise from the nonlinear response of 
the vacuum due to virtual electron-positron pairs. This 
can lead, for example, to laser photon merging [33| or 
the emission of Unruh radiation [39|, [40| . Further nonlin- 
ear vacuum effects in strong laser fields comprise photon- 
photon scattering [4l|, [43, \^ , photon splitting [4^, [4^ 
and changes in the refractive index [4y, [47|] . The creation 
of bound states in a supercritical ionic field was treated 
in HI. 

In the present paper, we consider the bound-free chan- 
nel of electron-positron pair production in XFEL- nucleus 
collisions. In contrast to the earlier treatments of this 
problem |36l . [37[, a laser beam of linear polarization is 
chosen. On the one hand this choice represents the more 
easily achievable case for experimental realization [l(l| 
but leads on the other hand to more difficult calculations 
due to the appearance of the generalized Bessel functions 
(see Eq. (|6])). Additionally, we consider not only the cap- 
ture of the electron into the ground state, but also into 
excited states in the L-shell. This gives rise to corrections 
to the total bound-free production channel of ~ 15 — 20%. 

The paper is organized as follows: Sec. [TTjdescribes the 
theoretical framework and the analytic calculations per- 
formed to obtain the result for the differential production 
rates. In Sec. IIIII numerical results for total and differ- 
ential pair production rates are shown and discussed. In 
Sec. IIVI a summary and conclusion are given. 

Relativistic units are used such that c = h = I and the 
elementary charge e = ^/a with the fine-structure con- 
stant a. The (scalar) four product is denoted, e.g., with 

(j)x) and p means p'^7^. The standard metric (H ) 

and Dirac matrices 7^ are employed [49j. 



Eq. ID) is the so called post form (in contrast to the 
prior form) of the transition amplitude. In the post form 
the bound state electron wave function is assumed to 
be free from the interaction for i — *■ cx3 while the free (i.e. 
unbound) positron wave function '5 feels both the laser 
and the Coulomb field. The laser four potential A'^ is 
considered as the interaction in this case and is turned 
off asymptotically at i — > 00. 

It should be noted that pair production is very similar 
to strong-field ionization [50, [Sl], [54] from a theoretical 
point of view. In ionization an electron is lifted from a 
bound state to a state in the positive-energy continuum 
via the absorption of photons. Similarly, in pair pro- 
duction an electron in the negative energy continuum is 
lifted into a bound state. Hence, initial and final states 
are essentially interchanged. While for pair production 
we describe the process in the post form, in ionization it 
is more appropriate to use the prior form [53l . [5J] . 

To be able to advance analytically one applies to 
Eg. (d]) the so called Strong- Field- Approximation (SEA, 
[55I . [Sg . [53, [5I, [lol), where one replaces the fully inter- 
acting positron wave function by a Volkov solution. The 
SFA can be applied when the laser field is sufficiently 
strong and the influence of the Coulomb potential on the 
positron is negligible. On the other hand the laser field 
strength must not be too strong so that the existence of 
bound states is assured. This sets limits on the nuclear 
charge number Z which should be of intermediate value 
to give reasonable results [36| . 

Within the SFA, the amplitude in Eq. ([T]) becomes 
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The so called Volkov states [60, l61| are solutions of the 
Dirac equation for a positron in a plane-wave field. They 
are defined as: 






"J^p.^explz^'+'l , (3) 



where p^ and k'^ are the four momenta of positron 
and laser photons, Vp,s are the conventional Dirac spinors 
(chosen like in |4£|) and the action reads 



5'-+^ = (px) + 
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p-A{fj)~-A^i^) df,. (4) 



II. THEORETICAL FRAMEWORK 

The process of bound- free pair production in combined 
laser and Coulomb fields can be described within an S- 
matrix formalism. The amplitude for the transition from 
the negative energy continuum to a bound state reads 

M 
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The normalization reads Np — \ P^ where q'' — p^ -\- 

22, 
-m§r^k is the effective four momentum in the presence 

of the laser field [6l| . 

In the present work, a linear laser polarization is cho- 
sen. The electric field points in the x-direction while the 
laser field propagates in the z-direction. Furthermore, 
temporal gauge is applied: 



Al = a(0, 1,0,0) cos ry 



A:^ =w (1,0, 0,1). 



Here, a is the amplitude of the laser field, to its fre- 
quency and r] = (kx) the laser phase. The action in 
Eq.(j4]) can now be written explicitly as 



/ 1 ^ ea e^c? 

5^+^ enters as a phase factor in Eq. ^. It can be 
expanded into a Fourier series, yielding the generalized 
Bessel functions of two arguments [53] : 



J„ (ai, 02) = ^ Jn-2m (ai) Jm ("2) (6) 



time evolution. C is a normalization constant, Z the 
nuclear charge number, as = 1/ (am) the Bohr radius 
and a — y/l — aZ the energy in units of the electron 
mass, i.e. Eis — ma is the (relativistic) bound energy 
for the Is state. The explicit forms of the spinor x* and 
C can be found in [3, [49| . 

In order to evaluate the pair production probability 
we follow the procedure outlined in [53l] for relativistic 
strong-field ionization of K-shell electrons (see also [3q). 
First we sum over the possible spin configurations of 
the electron and positron, thus assuming a polarization- 
insensitive measurement of the produced particles. Fur- 
thermore we have to calculate the square of the S-matrix 
amplitude according to 



The J„ denote the regular Bessel functions of first kind. 



e^a^ 



Their arguments ai = (f^yPa; and 02 = afi. \ 
prefactors in Eq. ([S]). 



are the 



A. Ground-state 

The calculation is performed in the ion rest frame, 
where the bound-state wave function in the ground state 
takes the form: 



'^is=5Mx^e-^^ 



Ei.t 



(7) 



It separates in a radial part g (r) — 
C(2Zr/aB)''"^e-^'"/''«, a spinor part x' and the 
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The sum over the electron spin is directly performed 
by adding the two spinor pairs M = ^^ Xs_Xs_ — 
X-1-1/2X+1/2 + X-i/2X'-.i/2 whereas for the positron we 
proceed by expressing the matrix elements of M in 
terms of 7 matrices, applying the well-known theorem 
J2s "^P^a = 2^ { p — to) . and calculating the traces via 
trace technology [43|. The time integral leads to the en- 
ergy conserving ^-function and the spatial integrals can 
be performed by elementary methods. 

The calculation naturally splits into three parts de- 
noted by the subscript letters A,B,C. The final result 
reads: 



J 
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UA + UB + uc] 



(9) 



Eg — q^ is the dressed positron energy, T — 2ttS (0) the interaction time [49| and no the minimum photon number 
due to energy conservation. The (5-function shows the energy conservation requirement explicitly. The term in the 
denominator represents the square of the Fourier transform of the (nonrelativistic) hydrogenlike wave function (j). It 
plays an important role in the spatial distribution of the emitted positrons (see Sec. IIIip . The functions ua,ub and 
Uc are defined as follows: 



UA = e^a^V (j,',_i 
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The arguments ai and a2 of the generahzed Bessel- 
functions J„ are suppressed here in order to maintain a 
better readabihty. L{,V,V are functions containing mo- 
menta, bound state energy and normahzation: 
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v = 



(l + a)(r(a)f22(-i) [l+(T) 



2-CT 



r(l + 2a) 



< OBpy 



Here, p = (f— nk = p — bk. Furthermore, r = ^—^ and 



~ 2(M ■ 

For the present process of pair creation, we are inter- 
ested in the production rate of electrons and positrons. 
Therefore, we have to integrate over the momentum 
of the free positron which has to be done numericaUy. 
It is advantageous to consider the effective momentum 

2 2 

qi^t — pM _l_ ^ ° fcM which includes the effect of the pon- 
deromotive motion of the positron. For small laser in- 
tensity the parameters p and q are approximately equal. 
The fully differential production rate reads 



d-'R 



T (2.f 



(11) 



Numerical results of Eq. pT|) are shown in Sec. IIIII 



B. 2si/2 and 2pi/2 states 

In this subsection the capture into the 2si/2- and 2pi/2- 
states is considered. The former gives the highest con- 
tribution to the pair production rate after the ground 
state, whereas the latter is the first state deviating from 
the spherical shape. 

The wave function for the 2s-state reads Q 



<I>L = 9 {r) X2s {r) e 



-iE2st 



(12) 



As before, g {r) — N2sr'^^^e^''^ contains the normaliza- 
tion and the radial part. Note that for the excited states 
also the spinor X2s depends on r. The new abbrevia- 
tion (^ = ''^^Xi^ ^^ used here. The bound energy reads 
E2s — 'm\/-^-^- The wave function of the 2pi/2-state 



'l±a_ 
2 

quantized along the beam axis is of the same form as 
the 2s-state with differences in the normalization and 
the spinor part [3| . The two states are degenerate in the 
Dirac theory, i.e. £'23 — E2p^/2 ■ 

The general procedure for the calculation of the 2s- and 
2pi/2-states is as before, but since the wave functions are 
more complicated, the final result for the production rate 
is more involved. We obtain 
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(13) 



The general form of the pair production rate (J13p is as in Eq. ([9]) . Once again, the energy conserving (5- function 
appears and the sum over all photon orders n is to be performed. Also, the square of the Fourier transform of the 
radial part of the Schrodinger wave functions appears in the denominator. The coefficients for the capture into the 
2s-state read: 



UA = e^a^P^^^) (^^ 
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and 



«C = ^^'c.2 (/ _ p^) p(2.) M ] ^2 M w^ + H2W2 + 2Hr7C-^ W3 (p. - bk^) \ 

X Jn-2^n-2 + 4J„_2 Jn + 2J„_2 Jn+2 + 4J„ J„ + 4J„J„+2 + Jn+2>/n+2 

Due to the energy degeneracy the final result, Eq. (fT3l) . holds also for the 2pi/2-state with the coefficients given by 
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New abbreviations used here are w = 



i + (p/cr 



ij = Till, H = Jig, 3i = 2i?, ^2 = -2!^ and la = ^^ = (^^^ + ^^C"' (^ + 1) rZ^)^ 



„ "j^. . The normalization factors are slightly different 



from before: 
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Furthermore, 



>v2-(gi + 2)v+32rVv)' 

W3 = pi + 2) 2iUV + pi + 2) J2C"^ (ff + 1) M^ 

+ lihC^UrV + llC^ {cr + 1) MrV 

Wa = {2iU + laC^ {<J + 1) mY 



We = Pi - 2) 3iZ^V + Pi - 2) n3C"^ (cr + 1) ,,WV 
+ lilsC^UrV + 2lC^ {a + 1) MrV . 

The abbreviations Z//,V and X are defined as before (see 



Eq. ([To])), additionally 



rU=il-(^] \wsmX + 2^wcosX 



and 

^V = cr (-2(cr + l)w + 3 + 2cr) sin A" 

- a ("2 (a + 1) ('^') I. - (a + 2) ('^') ") COS A- . 

Further analysis of the capture in the L-shell must be 
surveyed with numerical methods and will be discussed 
in the following section. 



III. RESULTS 

In the following, results for the numerical integration 
of Eqs. pT|) and ([13]) are shown. Laser parameters of the 
planned X-ray free electron laser (XFEL) at DESY are 
assumed (see [iO|)- If not otherwise stated, the intensity 
parameter is ^ = 10~^, frequency w = 9 keV, hence in- 
tensity / — 7.2 X lO^^W/cm^. Furthermore, the ions have 
a nuclear charge number oi Z = 50 and are accelerated 
to a relativistic 7-factor of 7 = 50. Primed quantities 
are in the ion rest frame whereas unprimed are in the 
laboratory frame. The total rate is not a Lorentz invari- 
ant since the relativistic time dilation has to be taken 
into account. The transformation law between the two 
frames simply reads 
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ab 



-Rir, 
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(15) 



In the multiphoton regime the process of lowest pho- 
ton order will dominate - for the parameters introduced 
above the minimal photon number is tiq = 2. Conse- 
quently, we will mainly discuss the two-photon process. 
The created electrons can be bound in the ground state 
but generally also in all higher states. Since the momen- 
tum spread is substantially larger for the ground state 
electrons, the capture probability to the K-shell is ex- 
pected to be the highest [3]. Capture into the L-shell 
yields the main correction to the total pair production 
rate. Hence, up to Sec. IIIIEI we will focus on the capture 
in the K-shell and work out the influence of the L-shell 
in Sec. ImFl 



A. One-Photon limit 

Pair production via the absorption of single 7-photons 
has been considered before (see, e.g., [3J,[3a|). The main 
goal of the present work is to regard nonlinear processes, 
thus exploiting the coherence of laser light. However, 
it is interesting to compare results obtained by the SFA 
methods of this work with earlier results of single photon 
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FIG. 1: Comparison of different calculations for the frequency 
dependence of the total pair production rate in the case of the 
absorption of a single photon (here, ^ — 10~^, Z = 1, 7 = 50). 



calculations, where the radiation field is treated pertur- 
batively. 

Figure [T] shows the scaling of the pair-production rate 
with the frequency of the laser for the absorption of a 
single high-energy photon (here, Z = 1). For this lin- 
ear process, i.e. the absorption of just one photon, sev- 
eral calculations are compared. A calculation in Born- 
approximation by Sauter |34l . l6l| and its high frequency 
limit are shown. We compare the frequency dependency 
of the total rate to SFA-calculations for linear and circu- 
lar polarization. One can see that the two polarizations 
yield the same rates as is expected for the absorption of 
a single photon. The main difference between the SFA 
calculation and the calculation in Born-approximation is 
that in the first case the positron is described by Volkov 
wave functions whereas in the latter the positron is a 
free wave. The two different calculation show a very 
similar frequency dependence and obtain the same high- 
frequency limit. Hence, the comparison of our calcula- 
tions to previously obtained results supports the applied 
methods and approximations. 



B. Scaling of total rate for 2-photon absorption 

In this section, results for the scaling of the total rate 
for the absorption of two photons are shown. 

Figure [5] shows the scaling of the total rate with the 
laser frequency for different atomic charge numbers Z. 
For tj = 6 — 16 keV the rate scales with {uj' — LUmin) 
for all Z, where (jJinin is the minimal frequency originat- 
ing from the energy conservation requirement. Above, 
the rate approaches a constant value due to the satura- 
tion of the available phase space. Below 6 keV, i.e. close 
to the energetic threshold, the rate's dependence on the 
frequency can still be described by a power law but the 
power is slightly larger than r = 1/2. It is interesting 
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FIG. 2: Dependence on frequency for various nuclear charge 
numbers Z. The hues are fits over calculated data points 
(here, i = 10"*, 7 = 50, n = 2). 
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FIG. 3: Dependence of the pair production rate on the atomic 
number Z for various frequencies. The lines represent fits over 
calculated data points (here, fioj = 6 keV (dotted), hio = % 
keV (dashed), Tiuj = 13 keV (solid), ^ = 10"'', 7 = 50). 



that the scaling-law exponent r depends on the laser po- 
larization. For a circularly polarized field the value r = 1 
was found [3a |. A similar behaviour has been derived 
for free pair creation by two-photon absorption [23|: In 
the case of linear (circular) polarization the cross section 
close to the energy threshold follows a power law with 
r = 2 (4). I.e., for both pair production channels, circu- 
lar laser polarisation leads to an exponent which is twice 
as large as for linear polarisation. 

Figure [3] shows the scaling of the total rate with the 
atomic charge number Z . A power law ^ Z'^ describes 
the calculated values well, we find a power d of 5.9, 5.7 
and 5.6 for photon energies of 6, 9 and 13 keV, respec- 
tively. For capture processes one generally expects a 
power close to 5 Q- The fitted curves in Figure [3] slightly 
deviate from the data points for small and high Z val- 
ues. When we consider only small Z, i.e. only weak 
Coulomb fields (hence small deviation of the SFA), the 
power approaches the value 5 in the high frequency limit. 
The same holds for the Born-approximation calculation 
by Sauter [34J|. Agger and S0rensen found in a calcula- 
tion for one photon pair production taking into account 
the full Coulomb effects a power slightly below 5 [3q . 

In summary, the final scaling equation for bound-free 
pair production by the absorption of 2 photons from a 
linearly polarized laser field reads: 



R oc 7-ie^Z5-^ {J 



.1/2 



(16) 



A similar scaling behavior exists for higher photon or- 
ders. The only difference is in the power of ^ - for the 
absorption of n photons it scales with ^^" . Since C <C 1 
in the multiphoton regime, the higher order terms are 
considerably smaller. 

For free-free pair creation the rate grows with Z^, 
hence the portion of bound-free pair creation becomes 
increasingly important for higher Z values. 



C. Angular distributions in ion frame 



The free positron is emitted under a solid angle dfJ' = 
d cos -d' dLp' . The distribution of the polar angle "d' is 
shown in Fig. [4] where d' = Q corresponds to the laser 
propagation direction. One can see that no positrons are 
emitted at i? = 0. The distribution has a peak at a small 
angle and falls off quickly thereafter. The general form 

n — ^ 

\2 



results from the term 



1 



in the equation 



{paB/zy 

for the final analytic result ([9|) . It represents the square of 
the Fourier transform of the (nonrelativistic) bound wave 
function. Compared are the distributions for higher pho- 
ton orders which have additional structure. The latter 
is a property of the generalized Bessel functions ^ and 
therefore a multiphoton effect. The higher photon order 
rates have been scaled by 10^^""^^ in order to be visible 
in this graph. The scaling of the rate with the intensity 
parameter is ^^" (see Sec. IIIIBp . 

For the case of a linearly polarized laser a dependence 
on the azimuth angle is expected (see Fig. [5]). The emis- 
sion is maximal in direction of the electric field compo- 
nent of the laser field and minimal but non-zero in direc- 
tion of the magnetic field component. For higher photon 
orders some additionally structure is present so that for 
odd photon numbers the minimum is shifted. The az- 
imuthal distribution is the same (up to total scaling by 
1/7) in the lab frame since the azimuth angle is not af- 
fected by the Lorentz transformation. 

The non-isotropic azimuthal dependence of the pair 
production rate is clearly a polarization effect and absent 
in the case of circular laser polarization considered before 
[36|. The positron distributions in the polar angle, how- 
ever, are similar for both polarization states. Regarding 
the total creation rates, we find that at the same intensity 
parameter ^ the circular case is favored, but when con- 
sidering the same intensity (/ — 7.2 x lO^^W/cm^), the 
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FIG. 4: Polar angular distribution of created positron in the 
nuclear rest frame (here, tko — 9 keV, ^ = 10""*, 7 = 50, 
Z = 50) for the absorption of two (solid line) , three (dashed) 
and four (dotted) photons. 



FIG. 6: Rate dependence on the energy in the laboratory 
frame for the bound-free (solid line) and free-free (dashed) 
channel (here, hco ^ 9 keV, ,f = 10"", 7 = 50, Z = 50, 
n = 2). 
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FIG. 5; Azimuth angular distribution of the created positron 
in the nuclear rest frame (here, fko = 9 keV, (, — 10~ , 7 = 50, 
Z = 50) for the absorption of two (solid line) , three (dashed) 
and four (dotted) photons. 



rate for linear polarization {R{i^ — 34.4 s^"'^) is slightly 
larger than for circular polarization (i?^;,.;. = 26.0 s~^). 
This comes from the fact that at same intensity the peak 
electric field of a linear polarized wave is higher (by a 
factor of \/2) . 



D. Transformation to the laboratory frame 

For experimental observation it is essential to consider 
the process in the laboratory frame. When one considers 
the polar angular rate or energy distribution the transfor- 
mation is no longer as simple as for the total rate. Instead 
Lorentz transformations have to be performed [3] and un- 
like before the energy of the emitted positrons is no longer 



fixed. In the ion frame, energy conservation fixes the 
positron energy to E' = 1.33 MeV. Now, for each angle 
i9' in the ion frame, there exists an energy Eg in the lab 
frame. Figure [6] depicts the energy-differential rate. It 
exists a minimal and a maximal energy corresponding to 
■d' — and tt and the distribution is peaked at Eq = 6.57 
MeV which matches the maximum of the fJ'-distribution. 
The positron energies are very high, the positron's rcla- 
tivistic 7- factor is 7+ — —^ > 10 so that the emitted 
positrons are highly relativistic particles. The distribu- 
tion for the free-free case is broader which is due to a 
broader angular distribution in the ion frame. 

Similarly, the polar angular distributions for bound- 
free and free-free pair production differ in the lab frame, 
as shown in Fig. [T] Basically all positrons are emitted 
under the same angle of z^min = 177.26° in the case of 
the bound-free channel. At smaller angles the emission 
of positrons is kinematically forbidden while at larger an- 
gles the emission is nearly completely suppressed. Picto- 
rially, the angles in Fig. [4]are squeezed into one due to the 
high velocity of the ion. This results from the form of the 
Jacobi-determinant that appears in the Lorentz transfor- 
mation. It is proportional to (l — /3//?+ cos??') , where 
(3 and P'_^_ are the relativistic velocities of ion (in the 
lab frame) and positron (in the ion frame), respectively. 
Since f3 > (3'_^, the determinant formally diverges at one 
angle zJmin, yielding the characteristic signature of the 
angular distribution. It should be noted however, that 
it is an integrable divergence not leading to difficulties 
in an experimental observation. The minimum angle 
"i^min coincides with the minimal accessible angle in the 
lab frame, sini^min — 7+/3+/7/3, which is described in 
[3|. For circular polarization the same angle t^min was 
found [33 • It is noteworthy that the angle of emission 
is unequal to 180° which would have been difficult for 
a possible experiment. Compared in Fig. [7] is the an- 
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FIG. 7: Rate dependence on the polar angle in the laboratory 
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gular distribution for the competing free-free channel. 
The positron distribution is also confined, but by far less 
than in the bound-free case. In an experiment this allows 
to distinguish both channels by positron detection only. 
This is important because the bound state can decay via 
photoionization. For the given parameters the two pair 
production processes have approximately the same to- 
tal rate (i?frco-frco = 0.98 S"\ ^bound-free = 0.69 S^^) 

but scale differently with Z and uj. To obtain an esti- 
mate of experimental pair production rates one should 
take into account the pulse length and repetition rate as 
well as the projectile beam density. For example, when 
A^ ~ lO^*' ions collide with a laser beam of T ~ 100 fs 
duration and / ~ 1 kHz repetition rate [lO|, then the 
number of pair creation events per second is of the order 
of RTNf/2 ~ ls~^ (assuming perfect beam overlap). 



E. High intensity domain 

For a relatively low intensity parameter ^ ^ 1, the pair 
production rate scales with f ^" . The lowest possible pho- 
ton order gives by far the largest contribution to the rate. 
If one increases the intensity, the higher photon orders 
become important whereas the two-photon process pe- 
ters out at some intensity (see Fig. |5|). The phase space 
for the two-photon process diminishes as one increases 
the intensity and at one point the process is no longer 
possible due to energy conservation: The ponderomotive 
energy of the positrons becomes too large. Eventually, 
the same happens successively for the higher photon or- 
ders but the black squares in Fig. [5] indicate that the 
sum over all photon orders continues to increase. To es- 
timate the contribution of all photons we calculated for 
each specific ^ the rate for photon numbers from 2 to 9 
and found out that they decrease exponentially for the 



larger n values. From the exponential fit we were able to 
perform the sum to infinity, hence giving us an approx- 
imate result for all photon orders. We point out that 
the effect described here is analogous to the phenomenon 
of channel closing in above-threshold ionization of atoms 

M- 

The results of this section are obtained by the same 
formalism as before; note, however, that the existence 
of bound states in the ^ '^ 1 domain is no longer evi- 
dent since the laser field strength is very large and even 
exceeds the atomic binding field. Therefore, the results 
should be considered with some care. 



F. L-shell contribution 

Apart from the capture in the ground state, the elec- 
trons can also be bound in higher states. Fig. [S] shows 
the distribution of the polar angle for different photon or- 
ders and the capture in the 2s-state. It is very similar to 
the distribution for the capture in the Is-state (Fig. [J]). 
The total rate, however, is substantially smaller. Table 
U shows absolute values of pair production rates for the 
Is- and 2s-states for different nuclear charge numbers Z. 
Their ratio 2s/ Is is nearly constant and approximately 
0.125. Furthermore, we find the same ratio for higher 
photon orders. Pratt [63] proved that for one-photon 
pair production in the high energy limit the ratio of the 
ns-shell to the K-shell is (T{ns) = a{K)/n^ (here, n de- 
notes the principal quantum number). Thus, we find that 
this also holds for higher photon numbers. 

The calculation of the 2pi/2-state is very similar to 
the 2s-state. However, the total pair production rates 
are substantially lower. Unlike before, for the capture 
into the 2p]^ /2-state no simple scaling law for the ratio 
to the ground state can be found. But whereas the 2s- 
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TABLE I: Total rates in the ion frame (in s ^) for different 
states and nuclear charge numbers Z (here, hu = % keV, 
i = 10-*, 7 = 50, n = 2). 
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FIG. 9: Polar angular distribution of the created positron for 
the capture in the 2s-state in the nuclear rest frame (here, 
fto; = 9 keV, C = lO"", 7 = 50, Z = 50) for the absorption of 
two (solid line), three (dashed) and four (dotted) photons. 



state scales with the same power of the nuclear charge 
number Z as the Is-state (^^ ■^^■^), this is no longer true 
for 2pi/2- Calculations for one-photon pair creation in 
Born-approximation suggest a power of Z^ and we find 
Z'^-^ . At Z = I the total rate of the 2pi/2 is a negligible 
fraction to the Is state (~ 10^^) but steadily increases to 
about 5% a,t Z — 92. Similar contributions were found 
by Agger et. al. for one-photon absorption [3^]; the 
2p3/2-state yields another ^ 2% at high nuclear charge 
numbers Z in this case. 

Fig. [TO] shows the (scaled) polar angular distribution 
for the capture into the Is-, 2s- and 2pi/2-states. One 
can see that the distributions are very much alike. At 
first sight this might appear surprising since we found 
in Sec. IIIICI that the general structure is governed by 
the terms in Eqs. ([9]) and (fT3|) which come from the 
square of the Fourier transform of the nonrelativistic 
bound state wave functions. When plotted as a func- 
tion of p' = (q'2 + n^Lj'^ - 2nJ(i cos d') ^'"^ the width of 
the Is state in momentum space is twice as large as the 
width of the 2s- and 2pi/2-states. But when plotted as a 
fmiction of i9', the width of the distributions turn out to 
coincide for all three states. 

In conclusion, the L-shell gives an additional contri- 



bution to the total bound-free pair production rate of 
15 — 20%. The next higher state, the 3s-state, can be 
estimated to increase the rate by 1/27 ~ 3%, according 
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IV. CONCLUSION 

In this paper we considered nonlinear bound-free 
electron-positron pair production in the collision of a 
relativistic bare ion with an intense X-ray laser field. 
The scaling of the total rate and the angular distribu- 
tions of the emitted positrons in the ion rest frame and 
the laboratory frame were analyzed. For the first time, 
bound-free pair production with a many-cycle, linearly 
polarized laser beam was investigated. Distinct polariza- 
tion effects were found in the azimuthal distribution and 
different frequency scaling. Furthermore, we analyzed 
the capture of the electron into higher atomic shells and 
found that these give a contribution of ~ 15 — 20% to 
the total bound-free pair production rate. The bound- 
free channel is comparable in the total rate to the free- 
free pair production channel for high ionic charges. In a 
possible experiment one should be able to distinguish the 
two channels by observing the angular distribution of the 
emitted positrons. For a laser like the planned XFEL at 
DESY the process could become observable when com- 
bining the laser with a sufhciently strong ion accelerator. 
Apart from the fundamental significance of the process, a 
possible future application would be to use the described 
setup for measuring high laser intensities via the total 
rates [3, 01 . Moreover, if anti-protons were used as pro- 
jectiles, the scheme might in principle be applicable for 
the creation of relativistic anti- hydrogen atoms [6J] . 
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